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The relativistic extension of the quantum harmonic oscillator, the so-called Dirac os- 
cillator (DO) [IH3], provides a paradigmatic and exactly solvable model of relativistic 
quantum mechanics. Originally proposed in quantum chromodynamics in connection 
to quark confinement models in mesons and baryons [3], the DO has received great 
interest in relativistic many-body theories and supersymmetric relativistic quantum 



mechanics (see [2j[3j|5HI] and references therein). The DO model is obtained from the 
free Dirac equation by the introduction of the external potential via a non-minimal 
coupling [ll[3l[7]. Since the resulting equation is linear in both momentum and posi- 
tion operators, in the nonrelativistic limit a Schrodinger equation with a quadratic 
potential is then obtained. In spite of the great amount of theoretical studies, the 
DO model in relativistic quantum mechanics and particle physics remains far from 
any experimental consideration. Recently, an increasing interest has been devoted to 
the investigation of quantum or classical systems capable of providing an accessible 
laboratory tool to test some peculiar phenomena rooted in the Dirac equation, such 
as Klein tunneling and Zitterbewegung [8]. Quantum and classical analogs of such two 
phenomena have been investigated in graphene, trapped ions, photonic crystals and 
metamaterials (see, e.g., [9HT5] and references therein), with the first experimental 
observations of Klein tunneling and Zitterbewegung reported in Refs. pMTT]. Such 
experiments motivate the search for experimental systems that could simulate a rel- 
ativistic DO. Recently, interesting insightful connections between the relativistic DO 
and the Jaynes-Cummings model of quantum optics have been pointed [ToT - fTS] . and 
an experimental proposal of a two-dimensional DO using a single trapped ion has 
been presented [TT] . 

It is the aim of this Letter to propose a classic wave optics analogue of the one- 
dimensional relativistic Dirac oscillator based on light propagation in engineered 
fiber Bragg gratings (FBGs). Let us consider light propagation in a FBG [19] with 
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a longitudinal refractive index n(z) = n + An h(z) cos[2tt z / A + (j){z)\, where n 
is the effective mode index in absence of the grating, An C no is a reference 
value of the index change of the grating, A is the nominal grating period defin- 
ing the Bragg frequency u>b = 7rc/(Ano), c is the speed of light in vacuum, and 
h(z), <j>(z) describe the amplitude and phase profiles, respectively, of the grating. To 
study Bragg scattering of counterpropagating waves at frequencies close to ujb, let 
E(z,t) = {u(z,t) exp(— zcugt + 27iizn / A ) + v(z,t) exp(— zcjgt — 2nizn / \ ) + c.c.} 
be the electric field in the fiber, where Ao = 2noA is the Bragg wavelength. The en- 
velopes u and v of counterpropagating waves [see Fig. 1(a)] satisfy the coupled-mode 
equations [T9|l20] 

[d z + (l/v g )d t ]u = iq(z)v (1) 
[d z -(l/v g )dt]v = -iq*{z)u (2) 

where q(z) = (irAn/ Xo)h(z) exp[i<f)(z)] is the complex scattering potential and v g ~ 
c/n is the group velocity at the Bragg frequency. To highlight the analogy between 
light propagation in a FBG with suitable amplitude and and phase profiles, and the 
one-dimensional relativistic DO, let us introduce the dimensionless variables x = z/Z 
and r = t/T , with characteristic spatial and time scales Z = Xo/(irAn) and T = 
Z/vg, and the new envelopes ^1,2(2) = l u ( z ) T v(z)]/y/2. In this way, Eqs.(l-2) can 
be cast in the Dirac form 

id T ip = a x {p x - if{x)a z } ip + a z m(x)ip (3) 
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for the spinor wave function ip = (ipi,if) 2 ) T , where a x and a z are the Pauli matrices, 



defined by 
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p x = —i(d/dx), and where we have set 



(4) 



m(x) = h{x) cos[0(x)] , f(x) = —h(x) sin[0(x)]. 



(5) 



In its present form, Eq.(3) is analogous to the one-dimensional Dirac equation [H], 
written in atomic units (K — c — 1), with a space dependence mass m and with the 
momentum operator p x substituted with p x — if(x)a z . The space dependence of the 
particle mass m is known to describe the particle interaction with a scalar Lorentz 
potential, whereas the substitution p x — > p x — if(x)a z corresponds to a non-minimal 
coupling which is essential to describe the relativistic DO (see, for instance, [7]). To 
realize the one- dimensional analog of the DO, let us choose the amplitude h and phase 
4> profiles of the grating such that hcos<f) = m and hsin<f) = —f(x) = u s m x, i.e. 



(see Fig.l' 



h(x) = m Q yl + (co s x) 2 , <j){x) = atan(w s x), 



(6) 



where m and u s are two arbitrary constants, corresponding to the particle rest mass 
and oscillation frequency of the DO in the non- relativistic limit [7]. Analytical expres- 
sions of the energy spectrum for the one-dimensional DO and of corresponding bound 



states can be derived following a standard procedure detailed e.g. in Ref. [7J. Let us 



search for a solution to Eq.(3) of the form ip(x,r) = (xjj + (x), ip-(x)) T exp(— iSr), 
where 8 is the energy eigenvalue of the Dirac Hamiltonian and m(x), f(x) are defined 
according to Eqs.(5) and (6). The functions ip + and ^- are then found to satisfy the 
harmonic oscillator equation 

1 d 2 i/j± 1 2 2 5 2 - ml =f m uj s 
- + 2^ ^ = 2^ ^ (?) 

The eigenvalues of Eq.(7) are those of the non-relativistic harmonic oscillator. As- 
suming for the sake of definiteness u s > 0, the positive-energy spectrum (electron 
branch) of the DO is thus given by 

S n = \Jml + 2m w s (l + n) , n = 0, 1, 2, 3, ... (8) 

whereas the negative-energy spectrum (positron branch) is given by 

5 n = —\J ml + 2m co s n , n — 0, 1, 2, 3, ... (9) 

The corresponding eigenfunctions ^i^) can be simply expressed in terms of 
Hermite polynomials multiplied by a Gaussian function. It is worth pointing out 
two important properties of the spectrum of the one- dimensional DO. First, the 
non-relativistic limit of the DO is attained by considering the positive-energy branch 
of the spectrum in the limit u s /m — > 0; in this limit, from Eq.(8) one obtains at 
leading order 5 n = mo + oj s (1/2 + n) + (u s /2) {n = 0, 1, 2, 3...). The first two terms 
appearing in this expression of the energy have a simple physical meaning: mo is the 
rest energy of the particle whereas oj s {1/2 + n) are the usual quantized energy levels 
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of the non-relativistic harmonic oscillator (please note that we used atomic units, 
i.e. h — c — 1). The third term (u s /2) is responsible for an additional energy shift of 
the zero-point energy, which is a characteristic feature of the DO (also found in the 
three-dimensional DO; see for instance [HIE]). A second property of the DO, which 
readily follows from an inspection of Eqs.(8) and (9), is that the negative (positron) 
energy spectrum is not obtained from the positive (electron) energy spectrum by 
sign reversal (5 —> —S), the positron branch possessing an additional bound state 
with energy 5 = —mo- In our FBG realization of the DO, bound-states with positive 
and negative energies should correspond to trapped light states in the FBG with 
resonance frequencies above (5 > 0) and below (5 < 0) the Bragg frequency u>b, 
respectively. In a FBG of finite length, the ideal amplitude profile h(x), defined by 
Eq.(6), must be truncated, i.e. one has h(x) = for \x\ > L/(2Z), where L is the 
grating length (see Fig.l). The effect of grating truncation is twofold. First, the 
bound states of the DO become actually resonance modes with a finite lifetime, 
which should be thus observable as narrow transmission peaks embedded in the 
stop band of the grating (see Fig.2 to be commented below). Second, the number of 
resonance modes sustained by the grating is finite owing to grating truncation. As an 
example, Fig. 2(a) shows a typical transmission spectrum [power transmission versus 
normalized frequency detuning 5 = (to — lub)T], as obtained by standard transfer 
matrix methods [19], of a FBG with length L/Z = 10 for parameter values m = 1 
and u s = 0.5. Figure 2(b) shows the corresponding reflection-band diagram (dashed 
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area) of the grating in the (x,5) plane [20], together with the electron (5 > 0) and 
positron (5 < 0) levels of the DO, as predicted by Eqs.(8) and (9); a snapshot of 
normalized intensity profiles \ip+\ 2 + \ tp-\ 2 = \u\ 2 + \ v\ 2 for a few low-order DO bound 
states are also depicted. According to Ref. [20], the band reflection diagram shows 
the local stop band of the grating, as a function of position x, which is determined 
by the inequality \5 — (l/2)(d<j)/dx)\ < \h(x)\. Note that the transmission peaks in 
Fig.2(a), embedded in the stop band of the FBG, occur precisely at the values 5 n 
predicted by Eqs.(8) and (9). Note also the asymmetry of the transmission spectrum 
around 5 = 0, with an additional resonance in the positron (5 < 0) branch of the 
spectrum with no counterpart in the electron (5 > 0) branch. According to the band 
reflection diagram of Fig.2(b), the number of bound states of the DO sustained by 
the truncated FBG are about 6 in each branch. For the sake of clearness, only 3 
modes in the electronic branch, and 4 modes in the positron branch, are depicted 
in the figure. Note that, as expected, the resonance widths increase for higher-order 
modes (i.e. for increasing quantum number n) because of the narrowing of the 
gap regions that sustain light trapping for such modes. Note also that additional 
oscillations in the transmission spectrum, outside the stop band of the grating, are 
clearly visible in Fig.2(a). However, such oscillations do not correspond to resonance 
states of the DO, rather they arise because of an impedance mismatch between 
the grating (\z\ < L/2) and non-grating (\z\ > L/2) regions of the fiber, like in an 
ordinary non-apodized FBG (see, for instance, [IE]). in physical units, for typical 
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parameter values n = 1.45, An = 1 x 10~ 4 and Ao = 1560 nm which apply to FBGs 
used in optical communications, the spatial (Z) and temporal (T) scales for the 
example in Fig. 2 are Z ~ 5 mm and T ~ 24 ps, respectively. Hence, in physical units 
the grating length is L ~ 5 cm, the unit scale of (non-angular) frequency detuning in 
Fig.2 is l/(27rT) ~ 6.6 GHz, and the maximum index change requested to realize the 
FBG, reached at the FBG edges z = ±L/2, is An^l + (co s L/2Z) 2 ~ 2.69 x 10" 4 . 
Such a nonuniform FBG should be easily manufactured with currently available FBG 
technology based on UV continuous laser writing [21]. According to Fig. 2(a) and 
following the previous discussion, in an experiment the resonant states of the DO 
could be simply detected from spectrally-resolved transmission measurements. Like 
in any FBG structure supporting trapped modes (e.g. Fabry-Perot or phase-shifted 
FBGs), the transmission peaks of high-Q resonant modes may be strongly reduced 
by unavoidable intrinsic losses in the FBG (see, for instance, [22]). As an example, 
Fig. 2(c) shows the transmission spectrum of the FBG of Fig. 2(a), but assuming a 
power loss rate a = 0.5 m _1 in the fiber (value taken from Ref. [22]). Note that, in 
spite of the great peak reduction, the resonant states as well as the asymmetry of 
the spectrum are still clearly visible. 

In conclusion, a photonic realization of a Dirac oscillator, based on light propagation 
in FBGs, has been proposed. Electron and positron bound states of the DO should 
be clearly visible in spectrally-resolved transmission measurements as resonances 
above and below the Bragg frequency. 
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List of Figure Captions. 

Fig.l. (a) Schematic of a FBG. (b) and (c) Example of FBG amplitude and 
phase profiles that realize the analog of the Dirac oscillator (m — 1, cu s — 0.5, 
L/Z = 10). 

Fig. 2. (Color online) (a) Numerically-computed power transmission (dB units) 
of a lossless FBG with amplitude and phase profiles shown in Fig.l, and (b) 
corresponding reflection band diagram (dashed area). In (b) a few low-order intensity 
profiles of trapped modes in the electron and positron branches are also depicted, 
(c) Same as (a), but for a lossy FBG (power loss rate a = 0.05 m^ 1 ). 
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Fig. 1. (a) Schematic of a FBG. (b) and (c) Example of FBG amplitude and 
phase profiles that realize the analog of the Dirac oscillator (m = 1, u s = 0.5, 
L/Z = 10). 
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Fig. 2. (Color online) (a) Numerically-computed power transmission (dB 
units) of a lossless FBG with amplitude and phase profiles shown in Fig.l, 
and (b) corresponding reflection band diagram (dashed area). In (b) a few 
low-order intensity profiles of trapped modes in the electron and positron 
branches are also depicted, (c) Same as (a), but for a lossy FBG (power loss 
rate a = 0.05 m _1 ). 
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